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ABSTRACT 

A software system is described which provides for batch least-squares estimation of spacecraft orbit, 
attitude, and camera bias parameters using image data from the Geostationary Operational Environ- 
mental Satellites (GOES). The image data are obtained by the Visible and Infrared Spin Scan 
Radiometer (VISSR) Atmospheric Sounder (VAS). The resulting estimated parameters are used for 
absolute image registration. Operating on the Digital Equipment Corporation (DEC) PDP-1 1/70 
computer, the FORTRAN system also includes the capabilities of image display and manipulations. 
An overview of the system is presented as well as some numerical results obtained from observations 
taken by the SMS-2 satellite over a 3-day interval in August 1975. 


217 



SECTION 1 - INTRODUCTION 


A variety of spacecraft (S/C) exist which transmit images to the ground to 
provide meteorological and Earth resource information. Several studies have 
been concerned with the use of this imaging data for the estimation of the S/C 
orbit and attitude. Such an estimation procedure can be used for several 
pui'poses. The one with which this report is concerned is the use of the esti- 
mated S/C orbit and attitude (O/A) parameters for absolute image registration. 
The estimated O/A parameters are used to predict the geodetic latitude and 
longitude (0, X) which correspond to a specified location on an Earth picture. 
This allows accurate geodetic coordinate determination for tempo ral phenomena, 
such as clouds or sea swells. 

There are two categories of image data; those from three axis stabilized S/C 
and those from spin-stablized S/C. 

The Landsat and Earth Resource Technology Satellites (ERTS) are examples 
of three axis stabilized S/C. These produce image data from high inclination 
(polar) close Earth (900 km altitude) orbits. The use of this data is discussed 
in Reference 1 which describes a software system for the display and manipu- 
lation of image data as well as the use of an extended Kalman filter estimator 
for the O/A parameter determination. 

The geosynchronous Geostationary Operational Environmental satellites 
(GOES) are examples of spin-stabilized S/C which produce image data. An 
overview of O/A estimation using this type of data is given in Reference 2, 
where sample numerical results are presented for the first geostationary 
Synchronous Meteorological Satellite (SMS-1). 



This paper describes a software system developed to provide Bayesian weighted 
least-squares estimation of spacecraft orbit and attitude parameters using 
picture data obtained from the VAS (VISSR Atmospheric Sounder) instrument 
to be flown on the GOES-D. The data consist of ground control points of 
known geodetic coordinates located on pictures of the Earth taken by the GOES 
spacecraft. The VAS/NAVPAK (VISSR Atmospheric Sounder Navigation Pack- 
age) system operates on the Digital Equipment Corporation PDP 11/70 computer. 
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SECTION 2 - VAS/NAVPAK SOFTWARE OVERVIEW 


As shown in Figure 1, the VAS/NAVPAK system can be divided into four 
functions. First, the Data Base Management (DBM) portion controls file 
and data manipulation. Second, the picture display and cursor navigation 
portion controls: (1) picture display on the I^S, such as image zooming; 

(2) cursor navigation, including the extraction of picture coordinates (X, e) 
and the amtomatic moving of the cursor to the picture coordinates correspond- 
ing to a specified longitude and latitude; (3) automatic grey scale correlation 
between a prestored chip (16 x 16 pixel reference landmark) and a search 
area about the cursor; (4) the creation of landmark observations. The third 
VAS/NAVPAK function is the o/A and camera bias estimation. This portion 
of the system provides for weighted least-squares PC) estimation of the 
satellite orbit, attitude, and camera biases. The fourth VAS/NAVPAK function 
produces the specific navigation parameters which are required over a spec- 
ified prediction interval (usually 2 days). The navigation parameters are 
used to annotate the picture data. 

2. 1 Picture Display and Cursor Navigation 

Cursor navigation is the prediction of picture coordinates (X, e) corresponding 
to a specified geodetic latitude and longitude, given the estimated satellite 
orbit and attitude and the camera biases for some epoch time. 

This is the method by which a prestored video reference area (taken from a 
VAS picture) is correlated with an area surrounding the cursor on the image 
displayed by the operator. 

2.2 OrMt/Atiitude Estimation 

The S/C O/A estimation is done with the classical Bayesian weighted least- 
squares technique. The estimator can use either landmark data, radar 
tracking data, or both. Only the capability for using landmark data will be 
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Figure 1. VAS/NAVPAK Overview 








presented. It is assumed that the working observation files of landmark 
data have been created before beginning the O/A estimation. 


The computational procedure followed for the O/A proceeds in the following 
steps: 


1. An a priori estimate is provided of the solve-for parameters. These 
parameters will be a subset of: 


Xj[> i=l, 5 
*i, i=1.5 
t 

s P 

at 0 


the S/C position and velocity 
S/C attitude model 
coefficients 
camera bias 
camera bias 
camera bias 


2. For each observation, the S/C position and velocity are found by 
integrating the equations of motion to the observation time, t Q ^ & . 

For the VAS/NAVPAK system, the integration is performed with a 
121h order Cowell method, as described in Reference 3. The force 
model is selectable by the user and can include a spherical harmonic 
geopotential expansion terms up to 21 x 21, lunar/solar third body 
perturbations, and solar radiation pressure. 

3. For each observation time, an observation (A ., e) pair and partial 
derviatives are computed corresponding to the geodetic coordinates 
{$, X) of the landmark using the S/C position, velocity, attitude, and 
camera biases. 


4. The computed observation pair is used to calculate the observation 
residuals. The residual is examined to see if it meets the editing 
criteria. If it does not, it is not used in the solution. 
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5. After steps 2, 3, and 4 have been performed for all the observations, 
the new estimate of the epoch S/C state, the attitude, and camera 
biases, and their covariance matrix, is computed. 

6, The new estimate of the solve-for parameters are compared with 
the previous to see if the least squares process has converged. If 
the solution is judged to have not converged, the new estimate re- 
places the a priori in step 1, and the process is repeated. 

2.3 Navigation Parameter Output 

Spacecraft parameters can be generated for a sequence of overlapping time 
intervals covering a specified output span. These parameters include 
spacecraft ephemerides, attitude information, camera biases, eclipse times, 
and Chebyshev coefficients for position, beta angle, and retransmission 
correction. 
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SECTION 3 - THE OBSERVATION MODEL 


The observational model in VAS/NAVPAK is a modification of that used in the 
SMS NAVPAK (Reference 4). The camera bias and attitude representations 
for the VAS/NAVPAK observational model were reformulated, consulting the 
VAS working group (Reference 5) and with the assistance of R. Pajerski (GSFC). 

The SMS and GOES are geosynchronous spinning spacecraft designed for talcing 
pictures of the Earth in several wavelengths. A camera, or VISSR (Visible 
and Infrared Spin Scan Radiometer), transmits data to a ground station where 
a complete picture of the Earth is assembled. The data consist of a grid or 
matrix of intensity measurements. A line number and an element number 
specify the location of the intensity measurement within the grid. The line 
number j 6, corresponds roughly to longitude. These are shown schematically 
in Figure 2. For the visible wavelength observations, each picture element 
(pixel) intensity measurement corresponds nominally to an area on Earth 
of dimension 1/2 mile by 1/2 mile square. Of course, near the edge of the 
Earth, foreshortening will enlai'ge and distort this square. Options exist to 
handle data whose dimensions are integer multiples of this unit (i.e. , 2-mile 
by 4-mile data). Associated with each line of the picture is a time and angular 
quantity which relates the starting position of the line to the direction of the 
Sun in inertial space. 

At the ground station preprocessing is performed and full resolution picture 
segments of 1024 x 1024 pixels are generated; In order to create a landmark 
observation, the operator first displays a picture or subset of a picture on 
the l^S. Then, an identification is made of a particular location on the 
picture e) pair which corresponds to a known geodetic latitude and longitude 
on Earth. The geodetic coordinates and the picture coordinates with associated 
quantities such as time and Sun angle are transferred to an observation file. 

This constitutes a single landmark observation pair. 
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Figure 3 shows the GOES satellite relative to the earth at an instant of time. 
Except for specific camera constants, the SMS is almost identical to the GOES. 
Both satellites are cylindrical spinning objects with the longitudinal symmetry 
axis nearly aligned with the spin axis. The spin axis in turn is nearly aligned 
with the polar axis of the Earth pointed southward. As the satellite spins, 
the camera scans across the face of the Earth's disk, from west to east 
measuring the light intensity for each pixel along a line. The relation between 
the (X, e) coordinates of each picture and the camera orientation can be shown 
by comparing the image in Figure 2 with Figure 3. The element, e, is re- 
lated to the azimuthal camera angle, q. This angle is measured in the satellite 
spin plane and is the angle between the line of sight (LOS) vector to the land- 
mark and the LOS vector to the left (west) edge of the Earth. The conversion 
to line element is 


e « q/RPE (1) 

where RPE is the number of radians per line element. The satellite spin plane 
in Figure 4, perpendicular to the spin axis z T , is shown coincident with the 
spacecraft (S/C) symmetry plane, perpendicular to the S/C longitudinal 
symmetry axis zg/ c . In the actual development of the observation equations 
the general case of a misaligned spin axis is considered. 


The line number related to the camera elevation angle, a, as 




RPL 


+ JLo 


( 2 ) 


where RPL is the number of radians per line and io is the line number which 
corresponds to a zero elevation setting of the camera. 


The relation of the picture coordinates (X, e) to coordinates of a location on the 
Earth (p, A) depends upon the spacecraft position and attitude, and the camera 
constants and biases. Several coordinate system transformations are required 
to express this relation. 
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Figure 3 


VAS/GOES Configuration 









The satellite spin plane coordinate system, the (x\ y', z') system in Figure 3 , 
must be related to the Earth inertial coordinate system in which (lie satellite 
position is computed. Figure 4 shows a spacecraft spin plane coordinate system 
relative to true-of-date coordinates. The x' axis lies in the true-of-date (xz) 
plane at an angle of x with respect to the true-of-date (-x) axis. The y' axis 
forms a right hand orthogonal system. 


The transformation matrix S from the (x, y, z) system into the (x', y', z') 
system is 


( 3 ) 


x -cosX 0 sinXl 

= S • | y = sin^sinX cosip cosXsiivp I 
z -cos^sinX sin ip -cos^cosXj 

Since the positive spin axis z' is nearly aligned with the negative z axis of the 
true-of-date system, the angles X and ip will always be relatively small. Also 
shown on Figure 4 are the right ascension and declination angles (a, 6) which 
are conventionally used to represent the location of the z' axis. The declination 
angle is near -90 degrees. The relation of (X, i]>) to (a, 6 ) is 


. v cosa 
tanX =^5 


( 4 ) 


sin^ =sinacos< 5 . 

The location of the z f axis in (a, 6) is expressed as a time varying function as 
6 — 6q + 6^t +62 sin(63t + 6^j) ( 5 ) 

and 


Q ' = + 0*1 1 + & 2 sin (<23 1 +cy 4) 


( 6 ) 


The model represented by equations ( 5 ) and (6) is a symmetric one. Because 
of the spin stability of the S/C axis, perturbations to (a, 6) or (x, ip) are ex- 
pected to be small. 
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NOTES: (1) (x, y, z) is the Iruc-of-dato system. The (xy) plane is the true- 
of-date Earth's equatorial plane* while z points along the Earth's 
axis of rotation. 

(2) (x\ y\ z') is the spacecraft spin system. The positive z 1 axis 
lies roughly in the direction of the -z truc-of-dnte axis. 

(3) The x f axis lies in the xz plane at an angle X above the -x axis. 
Hence* X is measured i ti the (x, z) (or, equivalently, (-x, z)) plane 
from the -x axis to the +x' axis. The angle X is measured posi- 
tive towards the +z axis. 


Figure 4. True-of-Date and Spin Plane Coordinates 
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Since the spin axis is nearly aligned with the negative z axis of the true-of- 
date system, the right ascension angle, a , is sensitive to the precision with 
which it is computed. For example, if the magnitude of the xy plane projection 
of z f (line OA in Figure 4) is nearly zero, then a change in sign would cause 
a to change by 180 degrees. Such a change can occur on successive iterations 
in the estimation process. The result would be to create divei'gent oscillations 
in the attitude correction vector (a, 6). Therefore, it is advantageous to use 
the (X, ip) coordinates for the spin axis location. 

The angle tp is analogous to declination and is the angle between the xz true- 
of-date plane and the spin axis. It is measured from the xz plane (perpendicular 
to the xy plane) to the z 9 axis. The angle X is analogous to right ascension 
and is the angle between the z axis and the projection of the spin axis onto the 
xz plane. 

The model for the (x, ip) coordinates of the spin axis can be written in a form 
similar to those of equations (5) and (6) 

# ~ + $4) (?) 


and 


^ X-jt + ^2 s ^ n (^ + ^4)* 


( 8 ) 


Figure 5 shows the S/C symmetry coordinate system. The z g ^ c axis is 
parallel to the S/C longitudinal symmetry axis and x g /£ points to the zero 
elevation angle in the actual VTSSR plane. The S/C symmetry plane is per- 
pendicular to the S/C symmetry aixs and is the reference plane from which the 
true camera elevation is measured. Two VISSR planes are shown; the actual 
VISSR plane is the plane swept out in elevation as the camera is moved from one 
spin cycle to another, while the nominal VISSR plane is the plane in which the 
camera motion is supposed to occur. The angle (measured in the symmetry 
plane) between the sun sensor plane and the nominal VISSR plane is y 0 . 
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(Xc/p, V q/ p, Z Q/r ) THE S/C SYMMETRY COORDINATE SYSTEM; Z S /C IS PARALLEL TO 
' 0/0 THE S/C LONGITUDINAL SYMMETRY AXIS, X S /c POINTS TO THE ZERO 

ELEVATION ANGLE IN THE ACTUAL VISSR PLANE. 

7 0 THE ANGLE IN THE SYMMETRY PLANE BETWEEN THE NOMINAL VISSR 

AND THE SUN SENSOR PLANE 

Am THE ANGLE IN THE SYMMETRY PLANE BETWEEN THE NOMINAL AND 

ACTUAL VISSR PLANES 

f IS THE ELEVATION BIAS ANGLE; THE ANGLE MEASURED IN THE 

ACTUAL VISSR PLANE BETWEEN THE CAMERA WHEN IT IS SET AT 
A ZERO ELEVATION SETTING AND THE TRUE ZERO ELEVATION 

a THE ELEVATION ANGLE OF THE CAMERA DURING THE nTH SPIN 

Aa THE ELEVATION INCREMENT DURING EACH SPIN. 


Figure 5. 


Spacecraft. Symmetry Coordinate System 





The angle measured in the S/C symmetry plane between the actual and nominal 
VISSR planes, ATo, is the first camera bias, A second camera bias is f, 
which is an elevation (or line) bias angle. When the camera is set at a zero 
elevation setting, represented by line x s / co , its true elevation angle is £. 

Angle a is the elevation of a landmark above (or below) the nominal zero eleva- 
tion point and Aa is the amount by which the elevation is incremented each spin 
cycle. At the beginning of a picture, a is set to a negative value corresponding 
to the northern part of the earth and then incremented to positive values towai'ds 
the southern portion of the earth. 

Figure 6 shows the spin plane (or attitude) coordinate system first instroduced 
in Figure 4. Because of the inertial motion of the spin axis (equations (7) and 
(8)) and the rotation of the earth, the location of a landmark with respect to the 
spin frame is changing. Moreover, the daily motion of die sun and the spin 
axis inertial motion causes the solar positions to change with respect to the 
spin coordinates. However, at the time of a landmark observation, t s , the 
azimuth of a landmark, and the azimuth of the sun, y^, can be determined 
with respect to the spin system. 

Figure 7 shows the spin coordinate system relative to the S/C symmetry frame. 
The symmetry frame is rotating with respect to the spin frame but Figure 7 
depicts the instant that the VISSR plane intersects with x 1 axis of the spin frame. 
Notice that the xg /q axis is shown coincident with the x f axis at this instant. 

This choice is tantamount to forcing the z f axis to lie in the Y^/q z ^/q pl an e. 

This choice as allowable because the bias f can absorb the elevation difference 
(between x r and xg/Q) which would occur if z' did not lie in the ys/c z S/C 
plane at this moment. 

An angle p is defined as the angle, measured in the spin plane, between the 
nominal VISSR and actual VISSR planes. This represents an azimuthal bias 
which allows the modeling of error in the azimuthal location of the VISSR plane. 


232 



/ 

I 

I 

t 

t 


X- 

<x\ y\; 
?2 


t 


2' SPIN AXIS 



:*) IS THE SPIN COORDINATE SYSTEM; 2' POINTS ALONG THE SPIN AXIS, Y' 
IS DEFINED SO THAT IT IS NORMAL TO THE SPIN AXIS AND LIES IN THE 
TRUE-OF-DATE XY PLANE, X' FORMS A RIGHT HANDED SYSTEM WITH 
Y' ANDZ\ 

THE AZIMUTH OF THE SUN IN THE SPIN P^ANE 

THE AZIMUTH OF THE LANDMARK IN THE SPIN PLANE. 

ISTHE ELEVATION ANGLE OF A LANDMARK ABOVE (OR BELOW) THE 
SPIN PLANE. 


Figure 6. Spin Plane Coordinate System 
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NOTES: 

Z' LIES IN THE Z $/c Y $/c PLANE 

a' IS THE ELEVATION ANGLE OF THE LANDMARK IN THE SPIN SYSTEM 

a + £ IS THE CORRESPONDING ELEVATION OF THE LANDMARK MEASURED IN THE 
ACTUAL VISSR PLANE 

£ IS THE CAMERA ELEVATION MISALIGNMENT BIAS. 


Figure 7. Camera and Spin Biases 
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The angle y is the angle in the spin plane between the actual VISSR plane and 
the censor plane, and a is the elevation of the camera at the time a landmark 
was observed in the S/C symmetry coordinate system. 

A A 

Since X'sun, Xsun and the sun direction form a right spherical triangle, 

4- - tan” 1 (tan A sin 6g) . ( 9 ) 

The picture coordinate, X, is the line number or elevation coordinate and is 
given by, 

1 = Kpt + £ o (10) 

where RPL is the radians/line conversion constant and is the line corre- 
sponding to a = 0. In practice the elevation angle a' is found in the spin plane 
and then converted to a. 

The second picture coordinate, e, corresponds to an azimuthal angle (measured 
in the spin plane) between the left (west) edge of the earth and the landmark. 

The situation is shown in Figure 8 which depicts the spin plane as viewed from 
the north. The satellite is spinning clockwise. The picture coordinate, e, is 
thus, 

7 2 - % - p + y -t 

e =— mod 2n (11) 

where RPE is the number of radians per element and P is the angle through which 
the satellite has turned from the instant of sun observation by the sun sensor 
to the observation of the left edge of the earth by the VISSR. The angle is 
the azimuth of the sun and Yg the azimuth of the landmark. The angle P is 
determined by finding, for each line of the picture, the first or leftmost pixel 
of that line. Each revolution, a body-mounted sun sensor on the satellite detects 
the sun and produces a sun pulse. For each revolution, a time interval called 
the (3-time (T^), is computed. This time, which should elapse between the 
sighting of the sun by the sun sensor and the alignment of the camera with the 
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EDGE 



THIS IS A VIEW OF THE SPIN PLANE SEEN FROM THE NORTH. THE (X' Y # Z' ) SPIN SYSTEM IS RIGHT-HANDED 
BUT APPEARS LEFT-HANDED IN THIS FIGURE BECAUSE THE SPIN VECTOR Z' IS POSITIVE INTO THE PAGE. 

7 1 IS THE AXIMUTH OF THE SUN {AT TIME t ) 

7 2 IS THE AZIMUTH OF THE LANDMARK (at t ) 

0 IS THE ANGLE THAT THE S/C HAS SPUN BETWEEN THE OBSERVATION OF THE SUN BY THE 
SUN SENSOR AND THE OBSERVATION OF THE LEFT EDGE OF THE EARTH BY THE USER. 

7 IS THE PROJECTION OF THE ANGLE BETWEEN THE VISSR AND SUN SENSOR PLANES. 

t s IS THE OBSERVATION TIME OF THE LANDMARK- 


Figure 8. Spin Plane 
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desired left edge of the earth picture, is used to detect the first element of each 
line. For each line, the values T g and t 0 (time of the average sun pulse) are 
available as recorded data. Since there are 3144960 counts per half spin 

' T ii 

3144960 (radianS) < 12 > 


or 


p 3144960 


n 


when the sun pulse is 180 degrees out of phase. 


(13) 
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SECTION 4 - NUMERICAL RESULTS 


The sample results shown below (Figure 9) are for a three day span of SMS-2 
data obtained from three images taken twenty-four hours apart. Additional 
preliminary results taken from NAVPAK runs using a longer data span supplied 
by NOAA indicates that sub-pixel accuracy is possible by using a suitable set of 
solve-for parameters and a longer, denser data set. The full results of these 
and other evaluations of VAS/NAVPAK (e.g. , force and attitude model evalua- 
tions, propagation/prediction capability evaluation, etc.) will be published in 
a future paper. 
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ITERATION REPORT FOR ITERATION 3 

CURRENT WEIGHTED RMS 0.727066D+01 PREDICTED WEIGHTED RMS 0.0000003+00 
PREVIOUS WEIGHTED RMS 0.7272O9D+01 SMALLEST WEICHTED RMS 0.727209D+01 

RELATIV CHANGE IN RMS 0.196389D-03 **** DC CONVRCED **** 


START = 750830 

1 50436 . 60 

END= 750901 150600 

.09 EPCCII= 750830 150000.00 



— 

OBSERVATION SUMMARY BY TYPE 


TYPE / 

TOTAL ] 

NO. / ACCEPTED / WEIGHTED RMS / MEAN RESIDUAL / 

' STANDARD DEV 

ELEM 

29 


24 

0 . 6566D+0 l 

-0.221 ID-01 

0 . 3283D+0 1 
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29 


24 

0.7913D+01 

-0. 1 195D-02 

0.3956D+01 

RANG 
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© , 0000D+00 

0.O000D+00 

0.00000+00 

RDIF 
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0 

0.O000D+00 

0.O000D+00 

0.00O0D+00 

ICEPLER I AN ELEMENTS AND 

LANDMARK MODEL ATTITUDE PARAMETERS 

FOR ITER 3 

PARAMETER 

SOLVE? 


CURRENT 

PREVIOUS 

STA.DEV 

SMA 

(KM) 



42164 . 9041 

42164.9041 


EGG 




0.0040 

0.G040 


INCL 

(DEG) 



1 .8162 

1 .8162 
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(DEG) 



128.2253 

128.2253 


CHI-1 

(DEG) 

YES 


0.4938 

0.4933 

0.3253D-01 

cm -2 

(D/S) 



0 f 0000 

0.0000 


PSI-l 

(DEG) 
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-1 .8674 

-1 .8674 

0.2964D-01 
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0.0000 
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CARTESIAN COORDINATES AND IARDMARK 

BIASES FOR ITER 

3 

PARAMETER 

SOLVE? 


CURRENT 
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STA . DEV 

X 

mi) 
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-26340.0797 -26340.0797 

0.2O0OD+01 

Y 

(KM) 
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32B82.0857 

32882.0857 

0.4616D+O1 

Z 

(KM) 

YES 


-778.1918 

-778.1918 

O.8142D+01 

XDOT 

(K/S) 

YES 


-2.4089 

-2.4089 

0.I779D-G2 

YDOT 

(K/S) 

YES 


-1 .9121 

-1.9121 

0 . 1 662D-02 

ZDOT 

(K/S) 

YES 


0.0790 

0.0790 

0.3193D-03 

BIAS-1 

(DEG) 

YES 


0.0970 

0.0970 

0.4514D-01 

BIAS-2 

(DEG) 



0.0000 

0.0000 


BIAS-3 

(DEG) 



0.0000 

0.0000 



Figure 9. Sample Numerical Results 
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